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The RHIC azimuth quadrupole: "perfect liquid" or gluonic radiation? 
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Large elliptic flow at RHIC seems to indicate that ideal hydrodynamics provides a good descrip- 
tion of Au-Au collisions, at least at the maximum RHIC energy. The medium formed has been 
interpreted as a nearly perfect (low-viscosity) liquid, and connections have been made to gravitation 
through string theory. Recently, claimed observations of large flow fluctuations comparable to par- 
ticipant eccentricity fluctuations seem to confirm the ideal hydro scenario. However, determination 
of the azimuth quadrupole with 2D angular autocorrelations, which accurately distinguish "flow" 
(quadrupole) from "nonflow" (minijets), contradicts conventional interpretations. Centrality trends 
may depend only on the initial parton geometry, and methods used to isolate flow fluctuations are 
sensitive instead mainly to minijet correlations. The results presented in this paper suggest that 
the azimuth quadrupole may be a manifestation of gluonic multipole radiation. 

PACS numbers: 13.66.Bc, 13.87.-a, 13.87.Fh, 12.38.Qk, 25.40.Ep, 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

The large elliptic flow v-i observed in RHIC heavy ion 
collisions compared to lower energies and to hydro pre- 
dictions is said to reveal a "perfect liquid" [l], Q • Elliptic 
flow is interpreted as a hydrodynamic response to early 
pressure and the azimuthal eccentricity e of the initial 
system Q. In the hydrodynamic (hydro) model of nu- 
clear collisions large flow values imply small viscosity and 
rapid thermalization Q. However, the applicability of a 
hydrodynamical model to nuclear collisions and the true 
nature of the azimuth quadrupole interpreted as elliptic 
flow can be questioned. In this paper I examine analysis 
techniques and interpretations leading to "perfect liquid" 
and describe some alternatives. 



A. Conventional flow argument 

The conventional flow argument is summarized as fol- 
lows 0, [B[ : Azimuthal asymmetry of participant nucle- 
ons at initial nuclear contact, copious particle produc- 
tion, subsequent thermalization by rescattering, and hy- 
drodynamic evolution lead to a matching asymmetry in 
final-state momentum space, providing direct evidence 
for strong interactions among initially-produced particles 
(gluons). Thermalization must be rapid because hydro 
expansion reduces the initial space asymmetry. The ar- 
gued connection between initial and final states through 
rapid thermalization is thus demonstrated by large ellip- 
tic flow and the apparent success of the hydro model, 
especially the relation vi oc e for ideal hydro. 

It follows that participant-geometry (eccentricity) fluc- 
tuations and flow fluctuations should manifest the same 
i)2 a £ relationship. Recent participant-based eccentric- 
ity models seem to unify flow measurements across colli- 
sion systems @ , and flow fluctuation measurements seem 
to provide an exact connection between flow measure V2 
and participant eccentricity e par t 0,9- However, model- 
ing low-x partons with a participant-nucleon eccentricity 



is questionable at mid-rapidity, and the true magnitude 
of V2 fluctuations remains unclear. 



B. Open issues for flow analysis 

Observation of the trend i>2 oc e would demonstrate 
that HI collisions at RHIC produce a thermalized QCD 
medium (quark-gluon plasma or variant) according to 
ideal hydro [H, [j] . Recent experimental results suggest 
that i)2fxe has been achieved [9(. However, there remain 
open issues: 1) Conventional flow measure V2, motivated 
by the hydro scenario 4], admits significant statistical 
bias (systematic error in numerical procedures), espe- 
cially for smaller event multiplicities, and is difficult to 
interpret [l(| ■ 2) Conventional flow measurements do not 
distinguish reliably between flow and "nonflow," identi- 
fied as angular correlations from low-Q 2 parton (gluon) 
fragmentation or minijets fill fl2L Ifli. 3) flow 
fluctuation measurements have been based on assump- 
tions about "nonflow" known to be invalid from minijet 
studies. 

4) Recent A-A eccentricity definitions are based on 
simulated distributions of point-like participant nucleons. 
Modeling the transverse distribution of low- a; partons 
at RHIC midrapidity with point-like nucleons is ques- 
tionable. Are assumptions and results consistent with 
known parton distributions, with the concept of a color- 
glass condensate? 5) Conventional flow measurements 
and terminology are strongly model dependent. Impo- 
sition of an a priori hydro model on data may obscure 
more fundamental processes. 



C. Possible resolutions 

Two new initiatives are needed. First, alternative 
physical mechanisms for azimuth multipoles should be 
explored. We should consider the analog in QCD to 
multipole radiation fields in electromagnetism, especially 
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given the non-Abelian nature of the gluonic field. What 
is the long- wavelength limit of QCD, the equivalent of 
Maxwell's equations? A static-field limit in the form of a 
color-glass condensate (CGC) has been introduced [ll|, 
but its dynamical implications have not been fully ex- 
plored. Gluonic multipolc fields may play a significant 
role in nuclear collisions and may produce strong azimuth 
correlations. 

Second, new correlation analysis methods include a 
model-independent statistical measure and analysis pro- 
cedure which reliably distinguish the azimuth quadrupole 
moment (model-independent terminology) attributed to 
elliptic flow from "nonflow" [1 01 ] . Accurate measurements 
of the azimuth quadrupole over a range of centralities, 
energies and nuclear sizes may provide new insights into 
the phenomenon described as elliptic flow. 

D. Paper outline 

In this paper I review conventional and 2D autocorre- 
lation analysis methods applied to measurement of ellip- 
tic flow or the azimuth quadrupole component, includ- 
ing power spectra, Pearson's normalized covariance and 
improved A-A centrality methods [lfj. I review two ap- 
proaches to A-A eccentricity simulation and advocate the 
optical Glauber model. I compare published STAR v 2 
data and fits to simulated 2D angular autocorrelations. 
I conclude that what is termed "nonflow" is dominated 
by the same-side minijet peak (jet cone) and that there 
may be a simple relationship between the true azimuth 
quadrupole moment and initial-state geometry parame- 
ters. I show that recently claimed large flow fluctuations 
are probably dominated by minijets. Finally, I review 
problematic issues for hydrodynamic models and suggest 
that the azimuth quadrupole is actually a manifestation 
of QCD field-field interactions, possibly related to longi- 
tudinal Alimentation and instabilities, that elliptic flow 
is actually gluonic multipole radiation. 

II. ANALYSIS METHODS 

Issues for azimuth correlation analysis include the cor- 
relation measure, the projection method to angular sub- 
spaces, the application of Fourier series, centrality mea- 
surement and plotting formats. Correlation analysis 
methods are also discussed in [lj], [III ■ 

A. Conventional method vs alternatives 

Conventional flow analysis is defined in terms of single- 
particle density p on azimuth angle <f> relative to reaction- 
plane angle \t r estimated by event plane angle Flow 
measure v 2 associated with the m = 2 Fourier term is ob- 
tained by several methods (e.g., event plane, subevents, 
two-particle correlations, four-particle cumulants) [Icl |. 



Because of a factorization in the definition of the Fourier 
series [l9j v 2 is actually the square root of the ratio 
of two pair densities, a per-pair two-particle correlation 
measure (lOj with the form v 2 = v/ Ap[2]/ p re f, where 
Ap = p — p re f is the density of correlated pairs, p re f is 
the reference (mixed) pair density, and '[2]' denotes the 
"second harmonic" or second Fourier coefficient of Ap. 
The result is a two-particle azimuth correlation analysis 
directed to a specific sinusoid amplitude. 

The event-plane method uses event- wise flow vector Q 2 
to estimate the reaction plane by event-plane angle ^ 2 . 
Intermediate value v 2 , observed 

= (cos(2[0 - \ff 2 ])) is then 
corrected by an estimate of the event-plane resolution. 
The EP procedure is closely related to a standard ID az- 
imuth autocorrelation analysis measuring the equivalent 
of v 2 {2}. Differences between v 2 {EP} and ^{2} arise 
because the EP procedure is an approximation to the ID 
autocorrelation [10| . 

The 2D quadrupole moment analysis is based on the 
power spectrum derived from the azimuth autocorrela- 
tion by a Fourier transform according to the Wiener- 
Khintchine theorem [lfj. When extended to 2D angu- 
lar autocorrelations on difference axes (77a,</>a) the au- 
tocorrelation analysis is able to distinguish "elliptic flow" 
accurately from so-called "non-flow" dominated by mini- 
jets [2(|. The four-particle cumulant «2{4} developed 
to suppress nonflow [21| approximates results from a 2D 
autocorrelation analysis. 

B. Correlation measures and Pearson's covariance 

The per-particle measure of azimuth correlations 
Ap[2]j \j pref has the form of Pearson's correlation co- 
efficient or normalized covariance converted to a den- 
sity. The '[2]' denotes the quadrupole component of az- 
imuth correlations. Its relation to v 2 is defined below 
and in [Iol |. The interpretability of per-particle corre- 
lation measures, their success in revealin g th e details of 
nuclear collision dynamics, is described in [12|, [H, [3 H2| ■ 

Pearson's covariance is invariant under linear super- 
position. Its value for a combination of independent and 
equivalent systems is the same as that for any component. 
Thus, changes in the measure under composition (e.g., A- 
A centrality variation) indicate real physical changes in 
the composite and truly differential correlation measure- 
ment. Because v 2 contains an extraneous factor 1/ \/n c h, 
variations with energy and centrality can be misleading 
(e.g., cf. 0). 

C. 2D angular autocorrelations 

Two-particle densities are defined on 6D momentum 
space {pt\,rii,4>i,pt 2 ,r]2,4>2)- In this paper I restrict to 
p t -integrated pair distributions on 4D angular subspace 
(171, 01,772) 4>2)- Further projection of the 4D space to 2D 
subspace (t]a,<Pa) (e-g- ?/a = r)i - m, = Vi + V2) 
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is described in [10|, LL2( ■ Within a limited rapidity inter- 
val near mid-rapidity stationarity is usually valid: the 
correlation structure is approximately uniform on pair 
mean pseudorapidity rjs/2, and we expect stationarity 
to hold on azimuth for a 2ir acceptance. Projection by 
averaging of the 4D space to difference axes (t7a,</>a) — 
a 2D (joint) angular autocorrelation — then discards no 
information, and correlations are undistorted by the pro- 
jection 

In contrast, projection to ID </>a or <j) — ^ r - 
reaction-plane angle), as in conventional flow analysis, 
abandons a large amount of information. Different dy- 
namical processes (e.g., flow and nonflow) confused in a 
ID azimuth projection may be easily distinguished with 
a 2D autocorrelation. Examples are shown in Sec. IIVI 



D. Fourier analysis 



defines an unbiased version of u|. Alternative plotting 
formats incorporating Ap\2]/ ^p re f and several geometry 
parameters are considered below. 



E. Azimuth quadrupole y z dependence 

The advantage of a per-particle measure is illustrated 
by the pseudorapidity dependence of «2 over a large n 
interval [24| • In Fig. [T] (left panel) pt-integrated data for 
Au-Au collisions at ^snn = 200 GeV are plotted in the 
conventional format t>2 vs r\. The shape corresponds to 
reported minimum-bias measurements (the solid curve 
sketches the data trend), and the peak amplitude agrees 
with the 40-50% [y ~ 4.3) centrality bin. The rectangles 
represent the STAR TPC acceptance. Flow appears to 
have a substantial magnitude even at large pseudorapid- 
ity, but errors there are typically consistent with zero. 



In conventional flow analysis a Fourier series represents 
ID event- wise azimuth density p((f) [n| [23|. In [lj| the 
mean value of the azimuth density po was factored from 
the Fourier series 
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p(4>) = Po { 1 + 2 } cos(m[4> - # r ]) 



(1) 



ty r being the reaction-plane angle. The quantities v m 
are then ratios of true Fourier coefficients. The ratio for- 
mulation makes physical interpretation difficult, in part 
because different dynamical processes contribute to nu- 
merator and denominator. 

v m seems to describe a single-particle distribution, but 
is not observable. Instead, event-plane (EP) angle 
^ m is estimated with particles from the collision, and 
the result is a two-particle correlation analysis which es- 
timates n(n — l)v^ n , nv^ or depending on algebraic 
details. Conventional event-plane flow measure 1>|{EP} 
approximates vfl^} from an exact sinusoid fit to the ID 
azimuth autocorrelation. 

In contrast, all dynamical processes can be represented 
by a 2D angular autocorrelation without invoking a phys- 
ical model. The per-particle angular autocorrelation on 
(t]a, 0a) is 03 



Ap 



Ap„f 



yPref \/Pre] 



(?7A,<M) + 2 2^ -=cos(m(M),(2) 



where the first term on the RHS describes "nonflow," 
that is, peaked structures on t]a and </>a, and the sec- 
ond term describes azimuth sinusoid components (mul- 
tipoles). The m = 3,4 terms may be measurable but 
are not essential to understand basic nuclear collision 
dynamics. We obtain a model-independent per-particle 
quadrupole component Ap[2]/ p re f rather than per-pair 
"elliptic flow" measure V2 motivated by specific hydro- 
dynamic expectations [H, . The relation between the 
measures 




Ap[2]/^fp~^] = n(n - l)(cos(20 A ))/27rn = n/2?r • ^,(3) 



FIG. 1: Left panel: Distribution of conventional per-pair flow 
measure V2 on pseudorapidity r\ for Au-Au collisions at 200 
GeV (trend adapted from 24]). Right panel: The same results 
converted to per-particle quadrupole measure Ap[2]/^/p re / 
and plotted on pion rapidity y z assuming (pt) ~ 0.45 GeV/c 
(solid curve). The dotted curve is a gaussian and the dashed 
curve is a beta distribution, the widths matching the solid 
curve. The dash-dot lines denote the STAR TPC acceptance. 

In Fig. Q] (right panel) v-i is converted to Ap[2]/y^p~^J 
(solid curve) using Eq. §S§ and the measured n — > dn/dn 
trend from [25| . The conversion rj — > y z assumes pi- 
ons with mean p t ~ 0.45 GeV/c to suggest kinematic 
limits. The quadrupole component falls to zero within 
the kinematic limits, and the statistical errors are ap- 
proximately uniform on rapidity, another advantage of 
a per-particle measure. The Ap[2]/y/p re f curve is well- 
described by a gaussian with a Vz — 2 (dotted curve), but 
the gaussian tails extend beyond the kinematic limits. A 
beta distribution consistent with those limits is plotted 
as the dashed curve [To] ], with parameters p = q = 4.8, 
r.m.s. width 1.7 and half-maximum points at ±2.4. 



F. A- A centrality measurement 

Significant improvements in A-A centrality determina- 
tion were introduced in (26| based on the approximate 
power-law form of the minimum-bias distributions on 
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rich (observed particle multiplicity in an acceptance) and 
on n part (participant number) and n^ in (binary-collision 
number). Accurate Glauber parameters n par t/2 and nu n 
are denned by running integrals on fractional cross sec- 
tion a /(To from N-N to b — A-A collisions. Accuracy 
for peripheral collisions is greatly improved by impos- 
ing extrapolation constraints (invoking information from 
measured p-p collisions). The combination provides cen- 
tralities accurate to better than 2% for all geometries. 

Centrality is measured by participant path length v = 
^ribin/npart, the mean number of N-N encounters per par- 
ticipant nucleon pair. Initial-state processes (e.g., parton 
scattering and fragmentation, minijets) measured per- 
participant should vary linearly with v for N-N linear 
superposition (l2l. [l3j. 



III. A-A ECCENTRICITY 

Interpretation of the A-A azimuth quadrupole requires 
an accurate eccentricity model, but there are major un- 
certainties about the correct model for nuclear collisions. 
The relevant model depends on parton x. Low-x trans- 
verse parton structure could be described by a contin- 
uum distribution, by point-like participant nucleons or 
by something intermediate. 

Eccentricity relative to the A-A impact parameter is 
defined as 



(y 2 ) - (* 2 ) 
(y 2 ) + (x 2 ) ' 



(4) 



where x and z define the A-A reaction plane [27] . The 
definition depends on the weighting function. In the op- 
tical Glauber model Woods-Saxon densities represent nu- 
clei A and B. The participant nucleon distribution is de- 
termined with thickness functions Ta, Tb- In a partici- 
pant Monte Carlo Glauber point-like participant nucle- 
ons are used for the weighting. In some cases the reaction 
plane is determined by the participant nucleons rather 
than the impact parameter of the parent nuclei [28l . |29| . 
The eccentricity then has two components 7r/4 apart 
(rank-2 tensors on azimuth) combined quadratically as 



{{y 2 )-{x 2 )} 2 + {±(xy)f 



^■part 



{(y 2 



2 >F 



(5) 



A. Optical Glauber eccentricity 

Some geometry parameters (6, n par ti v) are best es- 
timated by the participant-nucleon limit (Monte Carlo 
Glauber model) [26[. However, at small parton x the 
partonic azimuth correlation structure is the relevant is- 
sue for an eccentricity model. The optical Glauber model 
describes partons transversely correlated (bounded) only 
by the nuclear radius, a smooth distribution within that 
boundary. In the absence of other information that 



minimally-correlated configuration may be the best de- 
scription. 

In Fig. [2] the solid curves describe an optical Glauber 
estimate of 17 GeV Au-Au eccentricity with <jnn = 30 
mb described by the power series [3(| 

e(b) = -0.047 (b/b ) + 2.754 {b/b a f (6) 
- 4.797 (b/bof + 4.852 (b/b Q ) 4 - 2.492 (b/b Q f. 

Increasing (tnn to 40 mb causes a 5% reduction for 
central collisions increasing to 13% for peripheral colli- 
sions according to [3fJ, small compared to the dramatic 
difference between optical and participant Monte Carlo 
Glauber estimates. The hatched region in the right panel 
indicates the energy dependence of epsilon for peripheral 
collisions relative to the shape of the 30 mb estimate. 
We use the optical Glauber Eq. © (solid curves) in this 
paper to interpret quadrupole moments. 




FIG. 2: Left panel: Eccentricity e vs participant nucleon num- 
ber n pa rt modeled by an optical Glauber calculation (solid 
curve) and by a participant-nucleon Monte Carlo Glauber 
(dash-dot curve). The latter is a sketch of results in [tj [3l|. 
Right panel: The same curves plotted on mean participant 
path length v in number of encountered nucleons. The ver- 
tical dash-dot line denotes the mean value of v for N-N col- 
lisions [2(J . The hatched region indicates eccentricity uncer- 
tainties described in the text. 



B. Participant-nucleon eccentricity 

The dash-dot curves in Fig. [2] are obtained from a 
participant-nucleon (Monte Carlo) Glauber model and 
Eq. © [I, [H, [Hj]. The event-wise distribution of par- 
ticipant nucleons is used to estimate the reaction plane 
(which deviates randomly from the A-A impact parame- 
ter) [|| . The difference between optical and Monte Carlo 
curves is most dramatic for peripheral and central colli- 
sions. 

The large value e — > 1 for peripheral A-A — > N-N colli- 
sions implies that N-N collisions are on average rod-like 
(action at a distance). But nucleons are not point- like 
objects relative to the average interaction distance. Bar- 
ring contradictory evidence N-N collisions should be de- 
scribed by the average eccentricity (~ 0.25) of minimum- 
bias A-A collisions, since the geometry is in either case 
intersecting spheres. 
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Nonzero e for central A-A collisions implies structure 
resulting from modeling nuclei as distributions of point- 
like participant nucleons. What justifies that model for 
parton (gluon) interactions at x ~ O.Of? How does 
that model relate to the color-glass condensate (CGC) 
as a continuum limit? Since the number of partons 
is much greater than the number of nucleons at x ~ 
0.01 there may be substantial central-limit suppression 
of participant-nucleon correlation structure and fluctua- 
tions. 



IV. MEASURED AZIMUTH CORRELATIONS 

In [To| a simple model for p t -integrated flow centrality 
dependence was introduced. Flow and minijets (nonflow) 
inferred from 2D angular autocorrelations on ( 77a, 0a) 
were compared with ID projections on azimuth differ- 
ence variable 0a- The model was intended to illustrate 
qualitative features. In this section I pursue the compar- 
isons in more detail and examine novel plotting formats 
which reveal unexpected simplicity. 



B. Per-pair vs per-particle quadrupole measures 

Fig. [4] (left panel) shows per-pair V2 vs participant path 
length v, roughly proportional to the fraction of total 
cross section in the form 1 — a/ao (26j . The points are 
STAR measurements from [32j, and the dashed curves 



represent ID and 2D fits to the simulations described 
above. The hatched regions indicate the effect of the 
uncertainty in e for peripheral collisions shown in Fig. [5] 
(right panel). 

t>2 is substantial for peripheral A-A and N-N colli- 
sions and falls toward zero for central collisions, typical 
of per-pair correlation measures which contain an extra- 
neous factor l/n c h- The vertical dash-dot line marks the 
mean position on v of N-N collisions [lOj. t>2{2} "two- 
particle correlation" measurements arc typically larger 
than event-plane (EP or "standard" ) measurements [33| , 
although they estimate the same quantity. The i>2{4} are 
four-particle cumulant measurements intended to elimi- 
nate "nonflow" contributions (minijets) [2l[. The differ- 
ence between t>2{2} and t>2{4} has been interpreted as 
"nonflow" in the past, but more recently has been at- 
tributed entirely to vi fluctuations @, H[ • 



A. Fitting 2D angular autocorrelations 

The 2D autocorrelation model for this study is similar 
to [loj with two exceptions. Same-side peak amplitude 
0.6 v/2ir is larger by 15%, and the quadrupole amplitude 
model is Ap[2]/ ^/p^J — 0.0045 nun e 2 , as explained in 
Sec. ITVCl Fig. I shows examples for 90-100% and 20- 
30% centrality bins. 




FIG. 3: Left panel: Simulated 2D angular autocorrelation for 
200 GeV Au-Au collisions and 90-100% centrality (~N-N col- 
lisions). Right panel: autocorrelation for 20-30% centrality. 

Simulated 2D autocorrelations with statistical noise 
added (corresponding to about 5M events) for eleven cen- 
tralities were fitted with a model function consisting of 
constant offset, azimuth dipole, quadrupole and 2D gaus- 
sian same-side peak. 2D autocorrelations were also pro- 
jected onto 4>a and fitted with ID model cos(20a) to 
emulate a «2{2} analysis. Results for ID and 2D fits are 
discussed in the next subsection. 
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FIG. 4: Left panel: Conventional per-pair flow measures 
V2{2} and i>2{4}, with data from [§3] (points and solid curves) 
and corresponding simulations from this work (dashed curves) 
vs mean participant path length v. The vertical dash-dot line 
denotes N-N collisions [2a] . The hatched regions reflect un- 
certainties in the definition of e. Right panel: The same data 
and curves in terms of per-particle measure Ap[2]/^/p re /. 

In the difference between «2{2} and i>2{4} was 
comparable in shape and magnitude to the difference be- 
tween fits to the 2D angular autocorrelation and its ID 
projection (dashed curves). That difference is exactly 
the m = 2 component of the Fourier decomposition of 
the same-side minijet peak, implying that "nonflow" in 
conventional flow analysis is dominated by crosstalk be- 
tween minijet correlations and the quadrupole amplitude 
in a ID projection on </>a- 

Fig. [¥] (right panel) shows per-particle azimuth 
quadrupole amplitude Ap[2]/ y/p re f vs v. The dashed 
curves are ID (upper) and 2D (lower) fits to simulated 
autocorrelations, transformed from right to left panel via 
Eq. §5§ . The v 2 data from [32| (points) have been trans- 
formed from left to right also according to Eq. ([3]). The 
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quadrupole amplitude increases rapidly with centrality 
to a maximum for mid-central collisions, then follows the 
trend of the eccentricity toward zero. Although the cen- 
trality coverage of measured i?2{4} is limited, it is consis- 
tent with 2D fits to the simulated autocorrelations which 
extrapolate to zero at b = by construction. The agree- 
ment between simulations and data is good despite the 
simplicity of the model. 

This exercise was intended to demonstrate the effect of 
minijet correlations on v 2 {2}. The offset between dashed 
curves is due to the same-side 2D minijet peak, the sub- 
stantial m — 2 component in its ID azimuth Fourier 
expansion. The away-side minijet ridge does not con- 
tribute to V2{2} bias because it is described by dipolc 
term cos(</>a) and is therefore orthogonal to quadrupole 
cos(20a). The offset between i>2{2} and t>2{4} is con- 
sidered further in Sees. [V] and IVll An unanticipated 
benefit of the simulation is the agreement of measured 
Ap{A}/^p re f = rw|{4} /2ir with a simple n^ in e 2 trend, 
as described in the next subsection. 




n bin 1/S dn/dn 



FIG. 5: Left panel: Per-particle quadrupole measure 
Ap[2]/ y/p re f divided by initial-state geometry estimator e 2 
(optical Glauber) vs binary-collisions estimator nun (Monte 
Carlo Glauber), showing an approximately linear relation 
(dashed line) for Au-Au V2 {4} data at 200 GeV. Right panel: 
Conventional ratio U2/e vs in-medium collision-number esti- 
mator \/S drich/dr]. The hatched regions represent the low- 
density-limit (LDL) expectation (slope arbitrary) and the 
ideal hydro expectation. The dashed curve is the dashed line 
in the left panel properly transformed. 



C. Binary collisions vs hydrodynamics 

In Fig. [5] (left panel) the same V2 data are plotted in the 
format I/e 2 Ap[2]/^/p re / vs nbin, the number of binary 
collisions from a Monte Carlo Glauber simulation. The 
data derived from i>2{4} are well approximated by 

1/e 2 Ap[2]/^p^f = 0.0045 n bin , (7) 

suggesting that the azimuth quadrupole may be de- 
termined solely by initial-state collision parameters 
(b, ^/s~nn, A). However, accurate data are lacking, espe- 
cially for peripheral collisions down to N-N. The dashed 
line is the basis for the quadrupole amplitude in the auto- 
correlation simulations described in Sec. IIV Al The sur- 
prisingly simple linear relation suggests that the physi- 
cal mechanism of the quadrupole component is the same 
from N-N to central A- A collisions. It is therefore im- 
portant to test its validity with accurate data over the 
broadest possible centrality range and for other collision 
systems (e.g., lighter A- A and lower energies). 

In Fig. [5] (right panel) STAR v-i data are plotted in the 
conventional format V2/e vs 1/ S dn/dn, the latter reflect- 
ing the low-density limit (LDL) expectation that vije 
increases toward a thermal hydro limit with increasing 
number of in-medium particle collisions as part of an 
equilibration process [34j |. e from the optical Glauber 
model is used for this plot so that the two panels of Fig. [5] 
are consistent. Neither data set follows the hatched LDL 
trend (its slope is arbitrary). The dashed curve is trans- 
formed from the dashed line in the left panel. 

The i>2{4} data extrapolate to a value for N-N col- 
lisions which is more than 60% of central Au-Au colli- 
sions. Should we conclude that there is a strong tendency 
toward early thermalization and collective expansion in 
elementary hadronic collisions? We have no experimen- 
tal evidence for "saturation" at a hydro limiting value 



which might require an ideal hydro description. Prox- 
imity of data to a hydro prediction at one point does 
not imply that the hydro model is relevant to collision 
dynamics. Lack of extended agreement with either the- 
oretical conjecture contrasts with the simple relation in 
the left panel. 



V. NONFLOW VS FLOW 

Separating elliptic flow from other correlation sources 
(nonflow) is an unresolved issue for conventional f D flow 
analysis. Several strategies introduced to achieve sepa- 
ration [19t |21| rely on strong model assumptions about 
nonflow (e.g., hypothetical trends on A-A centrality and 
event multiplicity). Methods include higher-cumulant 
analysis (2f[ and flow- vector analysis [1, [n| [2l| . In (TTJ] 
we introduced azimuth multipole estimation from 2D an- 
gular autocorrelations which distinguishes different phe- 
nomena without physical model assumptions. 

A. Four-particle cumulants 

In [2l[ it was proposed that nonflow bias could be re- 
duced with higher-order cumulants. Higher cumulants 
(e.g., «2{4}) seem to reduce nonflow, but the model- 
dependent assertion that «2{4} eliminates nonflow con- 
tamination is questionable if nonflow is minijets as ob- 
served. The cumulant argument is based on the assump- 
tion that 1) nonflow is dominated by pair correlations, 
2) f;2{4} is insensitive to pair correlations and 3) true el- 
liptic flow is a universal characteristic of almost all (i.e., 
> 4) particles. Assumptions f ) and 3) are not generally 
true. 

Minijets are dominated at smaller p t by low-Q 2 partons 
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fragmenting to hadron pairs and triplets 20, 35]. But the 
assumption that minijets always contain less than four 
hadrons, especially at larger p t and more central colli- 
sions, is inconsistent with known fragmentation system- 
atica 35] and minijet measurements [12]. The assump- 
tion that elliptic flow is always "carried" by four or more 
particles is also unjustified, especially for more peripheral 
collisions where the quadrupole component may involve 
fewer than four particles. Thus, substantial bias may 
survive in t>2{4} measurements: positive bias (minijets) 
for more central collisions and larger p t , negative bias 
(number of "flowing" particles < 4) for more peripheral 
collisions. 



B. Variance of the flow- vector distribution 

Flow, nonflow and flow fluctuations have been stud- 
ied via systematic variations of the frequency distri- 
bution on "flow vector" magnitude Q 2 or normalized 
<l2 = Qil \fn (not to be confused with Fourier coefficient 
ratios q m = Q m /Qo defined in [l(|). Statistical reference 
Q le i is assumed to be gaussian distributed with variance 
. In the conventional description "nonflow" V n { 



"rcf " 

also contributes to Q 2 . Thus, the vector sum and Q 
are 



2 

ref 



h = Qref + V 2 +V n{ 



(8) 



Qref 



Ql - 2Q 2 V 2 cos(* 2 
2F 2 K f cos(* nf -* r ) + V; 2 f , 



where tildes denote event-wise random variables. To sim- 
plify I assume that flow V 2 is constant in magnitude, and 
nonflow is not correlated with the reaction plane (neither 
is true in general). The second sinusoid term is then zero 
in the mean, — > ng 2 in conventional notation, and the 
probability distribution on Q\ is 



dn 



1 



dQl ^h^; 



: exp 



V 2 2 + Q 



d 2 -^\, (9) 



a n = n (l + 52) being the variance of Q 2 about its mean 
V 2 , including nonflow 0, [lij [2l[ . The normalized vari- 
ance difference is ujl 



Acr^y = al/n - 1 = ng 2 /n = g 2 = nS 2 , 



(10) 



and Aa^/2ir = g 2 /2n (one unit of pseudorapidity) is the 
nonflow contribution to Ap[2]/ \/p r ef in Eq. ^j. 5 2 is the 
nonflow part of w|{2} = (cos(2[0a])) in conventional flow 
analysis. In a more general treatment V 2 may vary event- 
wise, and part of V^f may be correlated with the reaction 
plane. 

We have shown that 172 represents minijet correlations. 
"Nonflow" g 2 thus has a nontrivial dependence on cen- 
trality and multiplicity reflecting the physics of minijets, 
a complex subject undergoing intense study [20j. Various 



properties have been ascribed to g 2 . In [33[ g 2 is said to 
increase monotonically with centrality, contrary to previ- 
ous claims that nonflow is independent of centrality [13] • 
(In [33| the quantity g 2 — n part g 2 /n c h is plotted as g 2 , 
and in Eq. (10) g 2 — n part (w|{2} — i>|{4}) is defined.) 
However, no simple prescription can exclude minijet con- 
tributions from ID flow analysis. Conventional strategies 
result in an uncertain mixture of two physical phenom- 
ena. 

For this study I assume 1) "nonflow" cx minijets and 2) 
the centrality dependence is g 2 [v) oc v. The true minijet 
contribution to f|{2} is not strictly proportional to v\ 
it increases more rapidly because the same-side minijet 
peak broadens on ?7a with increasing A- A centrality [12| , 
thus increasing its ID projection onto 4>a relative to the 
same-side peak amplitude. I approximate the difference 
trend between dashed curves in Fig. [4] (right panel) by 



g 2 /2ir ~ 0.005 u, 



(11) 



with g 2 /2iT defined by Eq. (flU)) and the comment below 
it. 



VI. FLOW FLUCTUATIONS AND DATA 

Flow fluctuations have generated considerable recent 
interest. While the absolute relation of v 2 to hydro mod- 
eling may be uncertain, relative fluctuations could be ro- 
bust against such uncertainties and might confirm event- 
wise thermalization (v 2 oc e). By the same argument, if 
flow fluctuations are shown to be negligible then the con- 
ventional ideal-hydro flow scenario could be threatened. 
Thus, we examine flow fluctuations in some detail, ex- 
perimental aspects in this section and theoretical aspects 
in the next. 



A. Interpreting the f2{2} — W2{4} difference 

Recently, differences between v 2 methods previously 
attributed to nonflow have been ascribed instead to v 2 
fluctuations. From Fig. 2 of 0] we obtain 



v 2 2 {2} 
«2 2 {4} 
v 2 2 {2} - v 2 2 {A} 



a V2 and 



(12) 



v 2 — a 2 , therefore 



2a, 



where v 2 is nominally an unbiased mean value. How- 
ever, in Fig. 0] (right panel) we observe that the difference 
between A J o{4}/ v /p7^7 an d Ap{2}/ v /p7^J [n(u|{2} — 
v 2 {A])/2tt] corresponds almost exactly to the difference 
between Ap[2]/ \/p r ef obtained by a fit to a 2D angu- 
lar autocorrelation (lower dashed curve) and a ID fit to 
its projection onto 0a (upper dashed curve). Since that 
difference is clearly identified as minijet structure by its 
variation on 77A [101 ]. flow fluctuations inferred from the 
cumulant relations in Eq. (|12p are questionable. 
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B. Flow fluctuations from the Q\ distribution 

Eq. ((9]) is a conditional probability distribution func- 
tion or PDF: the distribution on Q\ given a fixed value 
of V 2 . If event- wise "flow" fluctuates then V 2 is also a 
random variable with its own distribution. The observed 
Q\ distribution is a folding of the RHS of Eq. ^ with 



dn 



1 



: exp 



2< 



(13) 



The PDF includes cr 2 ^ which broadens the PDF of Eq. © 
and which can in principle be obtained from fits to the 
PDF. However, nonflow also broadens the distribution on 
Q 2 ,, and any single fit to the distribution cannot distin- 
guish the two mechanisms. It is argued that the multi- 
plicity dependence of the two contributions at fixed cen- 
trality is different and can be used to separate flow fluc- 
tuations from nonflow. 

To examine the structure of the PDF on Q\ I simplify 
to ID, ignoring the out-of-plane component of Q 2 to fo- 
cus on gross systematic trends. Consistent with recent 
descriptions I rewrite the PDF in terms of q 2 = Qijypn. 
Ignoring constant factors the PDF on q 2 is given by the 
convolution integral 



the desired ideal-hydro relation v 2 oc e is obtained "by 
other means." Difference v 2 \2} — t>2{4} has also been 
attributed to v 2 fluctuations. However, the source of 
t>2{2} — W2{4} was identified by comparing 2D angular 
autocorrelations and their ID projections as in Fig. [4] 
which demonstrates that the v 2 difference interpreted as 
flow fluctuations is dominated by minijets. 

The apparent relation cr 2 2 /v| ~ cf/e 2 results from 
two misconceptions: 1) the trend identified as cr 2 2 /v 2 
reflects the relation in Fig. O (left panel) between the 
final-state momentum quadrupole component and the 
initial-state spatial quadrupole moment unrelated to fluc- 
tuations, and 2) the quantity identified as c 2 /e 2 models 
the small- a; parton distribution as point-like nucleons and 
NSD N-N collisions as point-like nucleons acting at a dis- 
tance. 

From Fig. [5] (left panel) we obtain Eq. ([7]) 
Ap[2]/ J~p^j ~ 0.005 nun e 2 , and from the t>2{2} — v 2 {4] 
difference in Fig. g] we obtain Eq. (fTT]) g 2 ~ 27r 0.005^. If 



we misidentify g 2 /2n 
factors we obtain 



a 2 " and ignore O(l) constant 



dn 
d% 



exp 



dv 2 exp 

(v 2 



fe ~ Vnv 2 ) 2 
I + .92W 



V2 



(14) 



2tt 0.005^ 
2n v\ 
0.005?/ 
2A (0 [2]/^7 
v 



(16) 



2< 



2 Tifjif 

1 



''part 1 



I let V2 vary about zero mean (v 2 — > 0) and obtain 



dn 
dql 



oc exp 



1 +92{v, n) + 2ncr 2 2 



(15) 



But since the participant-weighted eccentricity variance 
for Poisson statistics is "cr 2 " = cr 2 {n part } ~ l/n part it 
follows that 



ll_2 J! 



(17) 



It was assumed that any change in the width of v ^ e 2 

Eq. (|15p with random track discard isolates V2 variance 

cr 2 2 , based on the assumption that g 2 is approximately in- T h us , if the minijet contribution denoted g 2 /2n is inter- 



dependent of multiplicity [7j. However, systematic stud 
ies of minijet correlations show that two-particle correla- 
tions, particularly those represented by q 2 , decrease lin- 
early with random track discard [T2I l3q ] . Since the ex- 
pected width trend from random discard for v 2 fluctua- 
tions is 2n<7 2 2 and that measured for nonflow is <?2 oc n 
one cannot distinguish g 2 /2n from a 2 by random dis- 
card. "Flow fluctuations" inferred from the q 2 distribu- 
tion could as well be minijets, and minijet correlations 
may account for all such width variations. 



C. Flow fluctuations vs eccentricity fluctuations 

Recent studies of flow fluctuations emphasized the rela- 
tion of v 2 fluctuations and eccentricity fluctuations, and 
the trend cr 2 2 /w| ~ af/e 2 has been claimed 0, [1]. If 
v 2 fluctuations are equivalent to eccentricity fluctuations 



preted as "c 2 2 " and low- a: partons are modeled by point- 
like nucleons, v 2 data and Monte Carlo Glauber seem to 
support an event-wise v 2 oc e connection between initial 
and final states. 

In Fig. M (left panel "ct 2 2 "/u| - 0.015 vjnv\ (solid 
curve) results from confusing minijets with v 2 fluctua- 
tions. The r.m.s. trend is similar to that in Fig. 2 of 0]. 
The choice of V2 definition in the ratio (solid vs dashed 
curves) strongly affects the ratio distribution for periph- 
eral and central collisions. The hatched region sketches 
v 2 fluctuation data and errors from 0]. In Fig. [6] ( right 
panel) "er 2 "/e 2 ~ l/[n part e 2 } (solid curve) results from 
modeling transverse parton structure by point-like nu- 
cleons. Again, the choice among recent definitions of e 
strongly affects the r.m.s. ratio distribution for periph- 
eral and central collisions. For some choices the trends 
in the two panels are similar, giving the impression that 
V9 oc e |37| . 
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FIG. 6: Left panel: Relative r.m.s. 1)2 fluctuations vs A-A 
impact parameter b. The solid curve is derived from mini- 
jet trends and 1)2 obtained from 2D angular autocorrelations. 
The dashed curve is obtained using minijets and U2{2} from 
fits to ID projections. The hatched region represents the flow- 
vector analysis in Right panel: An estimate of the rela- 
tive r.m.s. eccentricity fluctuations from a participant-nucleon 
Glauber simulation (see text) divided by optical Glauber ec- 
centricities (solid curve) and participant-nucleon eccentricities 
(dashed curve). 



VII. FLOW FLUCTUATIONS AND THEORY 



sensitive to the Knudson number. 

The Knudson number Kn = n tot /n co u ~ X/ L is a mea- 
sure of internal collisions [38| . The reciprocal Knudson 
number Kn^ 1 measures the mean number of collisions 
per DoF during thermalization. Ideal hydro corresponds 
to Kn — > 0. In [38j a more complete vi fluctuation de- 
scription is proposed 



(18) 



Results in @, [|[ suggest that Jv\ ~ ^/ e " 
for mid-central collisions, and according to [38 
implying an upper bound on Kn ~ 0.03. 

That experimental upper limit is much smaller than 
UrQMD predictions, even after invoking artificially large 
cross sections to speed thermalization [38|. The results 
presented here suggest that the true upper limit on flow 
fluctuations may be much smaller than 0.15, and Kn may 
not be meaningful. 



According to some theory expectations flow fluctua- 
tions could result from fluctuations during system evo- 
lution toward equilibrium as well as from fluctuations 
in the initial geometry. Fluctuations could be gener- 
ated by in-medium collisions, the onset of turbulence or 
other aspects of thermalization. Fluctuation measure- 
ments might then provide access to equilibration dynam- 
ics and medium properties. 



A. Flow fluctuations and ideal hydro 

The main evidence for formation of a "perfect liquid" 
at RHIC is successful modeling of i>2 centrality and p t 
trends by ideal hydrodynamics [38] . In the ideal hy- 
dro model some initial conditions (e.g., spatial azimuth 
asymmetry) may be smoothly transported to a manifesta- 
tion in the final state, implying that fluctuations in the 
initial geometry could have a direct counterpart in the 
final state. Recent claims that relative V2 fluctuations 
are comparable to relative eccentricity fluctuations mod- 
eled by the participant nucleoli distribution 0, H[ seem 
to validate such expectations. 



B. Flow fluctuations and rescattering 

In a non-ideal hydro scenario other aspects of system 
evolution could contribute to flow fluctuations. Rescat- 
tering of a finite number of interacting quasiparticles with 
mean free path comparable to system size might con- 
tribute to «2 variance. If there is significant non-ideal 
system evolution (rescattering) V2 fluctuations might be 



VIII. GLUONIC RADIATION 

Given the initial and final states of a heavy ion colli- 
sion we seek the transport mechanisms which intervene. 
One source of information is propagation of the initial 
quadrupole moment to the final momentum space. The 
hydro model provides one description in the form of el- 
liptic flow. Its several parameters ought to have manifes- 
tations in the final state. Are they accessible in data, is 
the hydro model inevitable, are there alternatives? 



A. What can hydrodynamics describe? 

The hydrodynamic model applied to heavy ion colli- 
sions is summarized in 39]. The hydro sequence par- 
ton scattering — > fast thermalization — > flow with EoS 
— > hadronization provides one scenario for large-scale 
phase-space transport. However, hydro initial conditions 
are problematic: pQCD scattering cross sections are too 
small to thermalize initial-state partons [40j. There ap- 
pears to be insufficient time with known microscopic pro- 
cesses to unpack the nuclear wave function, equilibrate 
the results, flow them and reconstitute them into hadrons 
in the observed final state. In a viscous hydro model 
the viscosity-to-entropy ratio rj/s as a model parameter 
is driven to very small values in attempts to describe 
data ■ Does that imply a real medium with very small 
viscosity (perfect liquid) or an inappropriate model? 
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B. Implications from energy and centrality trends 

In Fig. [5] (left panel) 200 GeV Au-Au p t -integrated v 2 
data are consistent with 

(19) 

or n ch vl/2-n: w ttA v(b) {n part (b) e 2 {b)/2ir}, 

where the curly bracket on the RHS represent the per- 
particle quadrupole moment of the source, and the LHS 
is the per-particle quadrupole moment of final hadrons. 
Eq. (fl"9|) suggests that all p t -integrated quadrupole sys- 
tematics are described by one or two parameters rep- 
resenting the initial A-A system. There is no appar- 
ent sensitivity to intervening collision dynamics, no need 
to invoke a hydrodynamic scenario, equation of state or 
medium properties such as viscosity. The quadrupole 
may be completely determined by the initial small- a; par- 
ton (glue-glue) interaction. 

C. Implications from flow fluctuations 

Reliable separation of the quadrupole component from 
minijets reveals that flow fluctuations are much smaller 
than previously claimed. In Sec. IVICI cr^ 2 /i>§ was re- 
duced from mean value 0.15 pj to upper limit 0.03 for 
b ~ 7 fm. Fluctuations might be much smaller than 
expected from a participant Monte Carlo simply because 
there are many more gluons than participant nucleons — a 
central limit reduction. However, the absence of measur- 
able flow fluctuations may actually hint at the true trans- 
port mechanism, a simple relation between the hadronic 
quadrupole component and initial collision geometry de- 
fined by small- X gluonic field interactions. 

D. Generalizing parton scattering 

Parton interactions at large energy scales are modeled 
in pQCD as point-like interactions. Near the saturation 
scale, however, QCD interactions should extend over a fi- 
nite space-time volume — the "partonic participants" (in- 
teracting fields) may even extend across the nuclear di- 
ameter, the model implicit in [4ll |42|. 

In a generalization of the pQCD parton-parton vertex 
to non-pQCD interactions over extended space-time vol- 
umes the interaction strength should be the product of a 
cross section and a relative current density, a space-time 
current autocorrelation. Such an interaction would be 
consistent with Eq. I|19p provided an energy-dependent 
factor is incorporated. 

E. Gluonic multipole radiation 

Given present difficulties with the viscous hydro model, 
we can ask does a thermalized gluon population evolve 



according to the hydro scenario? Is a medium-dependent 
equation of state necessary, or even permitted by data 
trends? One possibility is that the hydro model is in- 
appropriate. The alternative may be an extended field- 
field interaction representing a generalization of pQCD. 
A radially expanding gluonic field might appear as a flow 
field, and some hydrodynamic properties (e.g., viscosity) 
could be ascribed to it. "Elliptic flow" might then be an 
instance of gluonic quadrupole radiation. 

In [42| and references therein interaction of chromo- 
electromagnetic (CEM) fields leading to Weibel (plasma) 
instabilities are studied in an attempt to solve the hydro 
initial conditions problem by replacing the pQCD 2^2 
partonic collision with extended field-field interactions as 
a route to fast local thermalization. But the CEM field 
approach could also supplant the hydro scenario entirely. 

The field-field interaction as a transport mechanism 
solves two major problems for quadrupole phenomenol- 
ogy: 1) the initial-final quadrupole relation Eq. ([7]) re- 
veals no information about the dynamical evolution of 
quasiparticles, and 2) flow fluctuations may be negligi- 
ble, implying either a large number of DoF or a contin- 
uum. Both problems are resolved if continuum field-field 
interactions provide the transport mechanism. The color 
filimcntation process may provide rapid large-scale trans- 
port from longitudinal to transverse phase space, and 
large azimuthal asymmetries may result [4lj . 

Consequences of an initial quadrupole moment have 
not yet been considered in studies of filimentation and 
instability. However, once CEM fields are the basis for 
dynamics multipole radiation is also a possibility. The 
final-state quadrupole in heavy ion collisions may be 
a manifestation of gluonic multipole radiation by anal- 
ogy with Maxwell's equations. Hadron azimuth corre- 
lations should reveal an approximately linear superposi- 
tions of gluonic multipole sources and corresponding radi- 
ated field components. The Weibel instability mechanism 
may in fact react to the initial participant eccentricity as 
a seed for instability growth, resulting in observed large 
azimuth asymmetries. 



IX. SUMMARY 

I have presented new analysis methods based on 2D 
angular autocorrelations with emphasis on azimuth cor- 
relations and the quadrupole component. I show that fits 
to 2D autocorrelations achieve model-independent sep- 
aration of "elliptic flow" (quadrupole component) from 
"nonflow" , and that the latter is dominated by the same- 
side minijet peak (jet cone). I have re-examine assump- 
tions about A-A eccentricity and conclude that estimates 
based on an optical Glauber model better describe low- a; 
parton correlations. 

Combining the optical Glauber eccentricity with pub- 
lished elliptic flow data I find a surprisingly simple linear 
relation to the number of binary collisions. If that rela- 
tion is confirmed by accurate data over a broader central- 
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ity range it is possible that "elliptic flow" is dependent 
only on initial-state collision parameters and is insensi- 
tive to hydro parameters or an equation of state. 

I also examine recent claims of large elliptic flow fluc- 
tuations comparable to participant Glauber eccentricity 
fluctuations and find that flow-fluctuation measurements 
are dominated by incorrectly-identified minijet correla- 
tions. True V2 fluctuations are much smaller, and present 
measurements are consistent with zero. 

Given those results and difficulties in accounting for 



hydrodynamic initial conditions (rapid partonic thermal- 
ization), combined with alternative theoretical scenarios 
invoking QCD field Alimentation and instabilities, I con- 
clude that the azimuth quadrupolc component may not 
be a hydrodynamic phenomenon, it may be an example 
of gluonic multipolc radiation. 
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